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Introduction 


The traditional stress tensor for a viscous fluid 

aui auj 
axj 

comprises terms representing the thermodynamic pressure, the 
volume viscosity, and the shear viscosity. Unlike the shear 
viscosity n, which can be derived from the Boltzmann Transport 
Equation and thus lends itself to a microscopic interpretation 
[1], the traditional volume viscosity \ has no physical basis 
and no relationship to established dissipative processes in 
fluids. In other words, neither its value nor even its order of 
magnitude can be predicted from the fundamental physical 
properties of a given fluid. In order to circumvent this 
obstacle to an orderly analysis of viscous flows, Stokes 
hypothesized that the volume losses in fluids are negligibly 
small [2] (which is true by definition for incompressible flows). 
This so-called "Stokes hypothesis" 

i = -(2/3) fi (2) 

has been utilized in a large volume of both experimental and 
theoretical data in many classes of compressible convective 
flows . 

In nonconvective (periodic or acoustic) flow, on the other 
hand, the point of view that volume losses are negligibly small 
has no basis in fact, to which the vast literature on sound 
absorption in fluids will attest. The underlying absorption 
processes in gases and some liquids are well understood and known 
to be attributable to a variety of relaxation processes, which 
accordingly constitute the physical basis for the volume 
viscosity. In gases the most prominent of these under ordinary 
experimental conditions are the relaxations of the molecular 
degrees of freedom. 

The purpose of this memorandum is to close the longstanding 
gap between acoustics and fluid dynamics with regard to 
volumetric losses in fluids. We start out by reviewing the 
fundamental physics of relaxation and its mathematical 
representation. Then we apply the lossy Navier-Stokes Equation 
to periodic (acoustic) flow and show that the traditional 
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expression for volume viscosity leads to a result which 
contradicts that describing acoustical relaxation. We 
demonstrate that the addition of a second volume viscosity term 
resolves the conflict and leads to a direct correspondence 
between the volume viscosity and the acoustic relaxation 
parameters. We proceed to discuss the formulation of volume 
viscosity in the presence of multiple relaxations. Finally, we 
present an example of the role of volume viscosity in a 
convective compressible flow. Although our discussion will be 
confined to gases, our analysis and conclusions will apply 
equally well to liquids since the mathematical representation of 
relaxation is the same in both media. 
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Re 
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real part of sound speed 

complex sound speed 

relaxed sound speed 

real part of acoustic wave number 

length scale for linearly accelerating flow 

length of contraction section 

Mach number 

thermodynamic pressure 

thermodynamic pressure at x=0 

relaxed pressure 

pressure increment due to kth relaxation 
reduced pressure, equation (37) 
amplitude of acoustical pressure variation 
Reynolds number, equation (32d) 
time 

time scale for pressure relaxation, equation (32c) 

one-dimensional flow velocity 

reduced flow velocity, equation (32a) 

ith component of flow velocity 

exit flow velocity of contraction section 

upstream flow velocity 

spatial coordinate 

ith component of spatial coordinate 
reduced distance, equation (32b) 
mole fraction of gas constituent 
imaginary part of acoustic wave number 
relaxation strength 
pressure relaxation coefficient 
density relaxation coefficient 
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K 

\ 

A 

P 

PW 

£° 

P(x) 

^R(X) 

A/> k ( x ) 

°ij 

T 

T ps 

T VS 


w 


reduced density relaxation coefficient, equation (32e) 

traditional volume viscosity coefficient 

acoustical wavelength 

shear viscosity coefficient 

gas density 

equilibrium gas density 

amplitude of acoustical density variation 
relaxed density 

density increment due to kth relaxation 
stress component 

experimental relaxation time, equation (lib) 
isentropic relaxation time at constant pressure 
isentropic relaxation time at constant volume 
(density) 

acoustical angular frequency 


Numerical Subscripts 


i,j,m indicates stress component (i,j,m = 1,2,3) 

k indicates relaxation process (k = 1,2,...) 

Mathematical Operators 


D/Dt 

= a/at + ua/ax, 


derivative 


= 0 if if*j, = l 

0( ) 

indicates order 


Acoustic Relaxation in Gases 


one-dimensional total time 

if i=j , Kroniker delta 
of magnitude 


The following derivation of the acoustic relaxation equations 
is based on small periodic (harmonic) variations, one-dimensional 
propagation, and a single relaxation process. The equations of 
motion 



£P 

ax 


and continuity 

au 1 dp 

ax p Q at 


lead to a quasi-wave equation 


(3) 

(4) 
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( 5 ) 


a 2 p _ a 2 p 
at 2 ax 2 


which is satisfied by travelling wave solutions of the form 
p = p e i (wt-kx) 

p — p e i(wt-kx) 


(6a) 

(6b) 


When we insert (6a) and (6b) into the acoustical equation of 
state [ 3 ] 


(1 + T ps st )0 - 77 (1 + T vs 

c o 


(7) 


*\j 

we obtain the ratio p/p, 
sound speed squared: 

£ 1_ _ J_ 

p c 2 c Q 2 1+ i WT ps 


which is the the reciprocal complex 

+ 0(t 2 ) (8) 


Upon separating equation (8) into real and imaginary parts, we 
isolate the terms representing the sound dispersion (real) and 
absorption (imaginary): 


c Q 2 _ w2r ps 2 

c 2 l+w 2 r ps 2 

or, equivalently 


(JT 


1 £ 


ps 


l+w 2 r 


+ 0 (£ 2 ) 


ps 


( 9 ) 


C 2 £ U 2 T 2 

= 1 + 

c 2 l+w 2 r 2 


+ i 


UT 




l+u z r 


2 ~ 2 


+ 0(f 2 ) 


( 10 ) 


As indicated in equations (9) and (10), a single relaxation 
process contains two defining parameters, a "relaxation strength" 
and a "relaxation time." The parameters used in phenomenological 
theories of relaxation are the isentropic relaxation times at 
constant pressure Tp S and at constant volume t vs , but those 
obtained from acoustical measurements are the relaxation strength 
£ and a third relaxation time r . The parameter pairs are 
interrelated by 
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(Ha) 


T ps” T vs 
€ = 

T ps 

T = J T ps T vs 

or conversely, 
Tps = T/J 1— € 

T VS = T J 1~« 


(lib) 


(12a) 

(12b) 


The relaxation strength is a measure of the capacity of the medium 
to convert translational energy of the flow into internal energy. 
It can be evaluated with great precision from known thermodynamic 
properties of the medium. The relaxation time is a measure of 
the rapidity with which the medium reacts to changes in 
translational energy. Because microscopic theories of the 
relaxation time are generally imprecise, quantitative evaluation 
depends on experiment. Information on the relaxation processes 
occurring in air are given later in this memorandum. 


It is important to note that the relaxation strength depends 
upon the difference between Tp S and r vs , a fact not taken into 
account in previous treatments of the volume viscosity. 

The real part of equation (10) yields the dispersion of the 
sound speed squared, shown in the plot of figure 1. The 
"dispersion step" has a height of f/(l-c), equal to the 
difference between the low and high frequency limiting solutions 
of the plot, and an inflection point where the condition ut= 1 is 
fulfilled. 


The imaginary part of equation (10) yields the sound 
absorption per unit wavelength: 


aA 



7T€ (JT 

J f 1+W 2 T 2 


(13a) 



(JT 


1+(J*T 


2-2 


aA 


+ 0 ( € 2 ) 


(13b) 



The plot of equation (13b) in Fig. 2 shows an absorption peak of 
height *€/(2 J l-< ) and a location at wr= 1. 

Both plots reveal that a relaxation process is ineffective at 
very low frequencies, where the internal degrees of freedom have 
ample time to equilibrate with translation, and at very high 
frequencies, where changes in translation occur too rapidly for 
equilibration to take place. In convective compressible flows, 
then, the volume viscosity would be expected to be most effective 
on time scales on the order of the relaxation time r . 


Traditional Volume Viscosity and Application to Acoustics 


We observe that the acoustical wave equation (5) does not 
contain a dissipative term. Dissipation is accounted for in the 
complex sound speed, which is derived from the acoustical 
equation of state (7) . The introduction of a physically 
meaningful volume viscosity implies a transfer of the relaxation 
terms from the equation of state to the Navier-Stokes Equation. 
Then, the equation of state will serve simply to determine the 
stationary temperature of the medium. 


To investigate the case of damped acoustical propagation, let 
us insert the traditional stress tensor [equation (1)] into the 
one-dimensional Navier-Stokes Equation, retaining only the 
dilatational terms. In place of equation (3) we obtain 


du _ do ll _ dp a 2 u 

p ° at ax ax + ,v ax 2 

where for one-dimensional propagation 
»; v = X+2 p 


(14) 


(15) 


Equation (14) together with the continuity equation (4) leads to 
the lossy wave equation: 


a 2 p *?V d 3 p d 2 p 

at 2 p o ax 2 at ax 2 


(16) 


Solutions of the form (6a) and (6b) yield the complex sound 
speed : 
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( 17 ) 


_ 1 _ k 2 

p c 2 w 2 -iw (»7 v //? 0 )k 2 

from which 

c Q 2 u> 2 ri v 2 /p 0 2 c 0 4 w>7 v //7 0 c 0 2 

= i + i + 0 (t ^ ) 

c 2 l+w 2 »7 v 2 /p 0 2 c 0 4 l+w 2 i7 v 2 /p 0 2 c 0 4 


(18) 


A term-by-term comparison between equations (18) and (9) reveals 
that the quantity t) v /p 0 c Q 2 can be identified with -Tp S , but that 
the two equations are compatible only if e=l — an absurd 
conclusion, for Kneser proves that the relaxation strength of a 
diatomic gas is bounded by e <0.082 for the vibrational relaxation 
and €<0.16 for the rotational relaxation [4]. Therefore, we 
conclude that the volume viscosity cannot be represented by the 
traditional constitutive equation (1) , containing only a single 
term. 


New Constitutive Equation for the Volume Viscosity 

Taking a cue from the acoustical equation of state (7) , we 
hypothesize that the second volume viscosity term is proportional 
to the time derivative of the pressure. Then the stress tensor 
becomes 



Dp au m 


aui auj 2 au m ' 

Q 

W- 

V_J. 

II 

1 

[ p - 'P Dt - a* m J 

6 ij + M 

axj axi 3 ax,, 


where for one-dimensional flow the factor 2/3 is replaced by 2, 
and the operator 


d _ a a 

Dt at ax 


is used to allow for pressure changes in convective flow. Let us 
examine the acoustical consequences. When we substitute the 
nonconvective dilatational terms of (19) into the one-dimensional 
Navier-Stokes Equation, we now obtain 



da ll 

ax 


ap a 2 u a 2 p 

ax + dx 2 + axat 


( 20 ) 


and the lossy wave equation becomes 
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( 21 ) 


d 2 p d 3 p _ d 2 p 3 3 p 

at 2 Po ax 2 at ax 2 p ax 2 at 


As before, solutions of the form (6a) and (6b) 
complex sound speed 

£ _ l _ k 2 -iwk 2 i? p x l-iwijp 

p c 2 u 2 -iw (i? v /p 0 )k 2 c Q 2 l-iw!j v //? 0 c 0 2 


lead to the 
+ 0 ( € 2 ) 


( 22 ) 


Comparison of the complex sound speed derived from the volume 
viscosity [equation (22) ] with that derived from the acoustical 
equation of state [equation (8) ] yields the volume viscosity 
parameters in terms of the acoustical relaxation parameters: 


*/p — ” T vs (23a) 

>»v = _ / , o c o 2 ' r ps (23b) 

Both »jp and r) v are negative quantities in order to provide the 
proper phase relationships when equation (21) is applied to 
acoustic flow. 


Volume viscosity and pressure relaxation cannot be modeled 
adequately via a single constitutive constant. The difficulty 
is due to the more fundamental character of pressure, in 
comparison to either volumetric viscous effects or shear 
stresses. Since pressure is accepted universally as a 
thermodynamic variable which can be used in the determination of 
other thermodynamic properties, it is not logical to force 
pressure to be related directly to the material rate of change 
of density, through the conservation of mass equation and a 
single volume viscosity coefficient. Furthermore, it is not 
logical to assume that pressure and density are determined 
uniquely for unsteady flows of simple fluids through an 
equilibrium equation of state, if those fluids possess internal 
degrees of freedom (which equilibrate at different rates) . These 
relaxation effects are easily visualized in terms of the behavior 
of air when it is subjected to high frequency acoustic 
excitations. It is certainly evident that the temporal behavior 
of the density, produced by these acoustic oscillations, does 
not correspond to that produced by a reversible, quasi-static, 
cyclic equilibrium process. In fact, it is known that pressure 
must respond more rapidly to these propagating disturbances than 
density because the disturbances consist primarily of collisional 
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exchanges of translational momentum, and pressure is a direct 
measure of translational momentum. Furthermore, since the 
nitrogen and oxygen molecules in air possess internal degrees of 
freedom that equilibrate more slowly than their translational 
counterparts, acoustically excited air is forced to exist in a 
regime where density and pressure are no longer in phase with 
each other, and temporal variations in internal molecular 
temperatures can be substantially different from their 
translational temperature histories. While it may be possible to 
ignore volume viscosity contributions completely, via Stokes' 
hypothesis, it is not possible to ignore these phase shifts 
between density, temperature (s) , and pressure that are manifest 
by these fundamental departures from thermodynamic equilibrium. 

If pressure equilibrium is enforced by an idealized equation of 
state, then the pressure is modified improperly via the 
traditional volume viscosity. It is certainly true that the 
pressure relaxation constant, introduced in our constitutive 
model, can be used along with volume viscosity; but, even 
further, it is quite possible that the pressure relaxation 
constant follows only the dynamic viscosity in importance as a 
constitutive constant. 

Multiple Relaxations 

The volumetric dissipation in many gases, notably air, is 
characterized not by one relaxation but by multiple relaxations 
occurring simultaneously. For the acoustical case Bauer shows 
that the effects of the individual relaxation processes upon the 
compressibility, as evidenced by the dispersion and absorption 
terms in equation (9), are additive [3]. The treatment here for 
generalized flow, then, is similarly based on the assumption that 
the effects of the individual relaxation processes are additive. 

A discussion of coupling among the relaxation processes appears 
at the end of this section. 

We define p R (x) and p R (x) as the "relaxed" pressure and 
density, i.e. the solutions obtained without losses. Then we 
assume that the pressure p(x) and density p(x) are composed of 
the relaxed contribution plus the sum of the contributions from 
all the relaxation processes. If the gas is a mixture, then each 
contribution must be multiplied by the mole fraction of the 
relaxing component of the mixture: 


9 



P(X) = Pr(x) 

+ l X k Ap k (x) , 
k 

k = 1, 2, ... 

(24a) 

p(x) = P R (X) 

+ l X k A/> k (x) 


(24b) 


k 


where Ap k (x) and Ap k (x) are the changes in pressure and density 
due to the kth relaxation alone. These changes are not 
necessarily perturbations and can, in fact, be very large. If 
Pj c (x) and p k (x) are solutions to Navier-Stokes Equation for the 
kth single relaxation process alone, using equation (19) , then 
the changes are simply 

APfc(x) = p k (x) - p R (x) (25a) 

A^k( x ) = Pk( x ) ~ Pr( x ) (25b) 

The above procedure is based on the condition that the 
presence of one relaxation process does not influence the effects 
due to another; in other words, the relaxations are decoupled. 

In reality this is not the case, for there are two types of 
coupling among the individual relaxations: gas-kinetic and 

nonlinear. 

Gas-kinetic coupling arises from the fact that the reactants 
and products of molecular reactions can appear in more than one 
reaction. The relaxation times Tp Sk and r vsk are no longer those 
of the decoupled kth relaxation alone. Rather they are 
determined from the solutions to a set of linear simultaneous 
equations involving all the coupled reactions. Since the 
coupling is linear, the formulation expressed by equations 
(24) -(25) remains valid; only the values of Tp Sk and r vsk have to 
be adjusted. For air this subject is treated in reference 5. 

The second type of coupling is due to the nonlinear nature of 
the convective operator u3/dx. When operating on equation (24) , 
it will introduce mixed terms in the Navier-Stokes Equation. 
Additional terms will be required in the series to account for 
the mixing. Generally, it may be expected that the volumetric 
losses and pressure and density changes are sufficiently small to 
permit equations (24) -(25) to remain very accurate approximations 
for most gases. 
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An Example: Linearly Accelerating Flow 


In order to examine the behavior of this constitutive model, 
it is convenient to examine a class of steady, one-dimensional 
flows. That class of flows includes flows through standing 
normal shock waves, but we are interested here in flows which are 
more typically encountered in the subsonic contraction sections 
of wind tunnels. If the velocity is given by u = u(x), the 
material or total time derivative is given by 

D d 

Dt dx 

and conservation of mass requires that 
d 

— (pu) = 0 (26) 

dx 

Thus, if we assume that the density and velocity are prescribed 
at the coordinate origin — say p G and U Q — conservation of mass 
is equivalent to the requirement 

P (x) u (x) = p 0 U 0 (27) 

Conservation of linear momentum can be written 


du 

pu— 

dx 


d_ 

dx 


(Po u o u ) 


do ll 

dx 


(28) 


and if we assume that 0 ^( 0 ) = -p D , equation (28) can be 
integrated to yield 

^il(x) = -p 0 - p 0 U 0 [U 0 -u(x) ] (29) 

which is a kinematical relationship not affected directly by the 
constitutive model. 

The constitutive model can be examined via the kinematical 
equation for normal stress. That is, must satisfy the 

relationship 


<>ll(x) = -p(x) + t ) pU — + r/ v — - -p Q “ p 0 U 0 (U 0 -u) (30) 

Before proceeding to the particular flow case, it is useful to 
distinguish between ''lossless' 1 pressure variations (where r/p, >?v' 
and p are taken to be zero) and "lossy" effects. That is, we 
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define Ap(x) : 


P(x) - p 0 - P o u o( u o" u ) 

Ap(x) (31) 

^o u o 

and introduce the reduced and dimensionless variables 


u = U/U Q 

(32a) 

X = X/L 

(32b) 

T — — ^p^q/L = Mc 0 t vs /L 

(32c) 

R e “ PoVqWp = P oM c o L //i 

(32d) 

K = -Vv/t* = p o c o 2r ps/^ 

(32e) 

i, our constitutive model must satisfy 
AP = -Tu (dAp/dx) + (Tu - /c/R e ) (du/dx) 

(33) 


For simplicity, we consider linearly accelerating flow 
systems, where: 

u(x) = 1 + x ( 34 

We note that in the case of this "contoured" wind tunnel 
contraction flow, the characteristic length L is related to the 
length of the contraction section L^ax anc * exit velocity U e 
by: 


L = 


kmax u o 

U e - u Q 


with 0 £ x £ 



(35) 


Hfire the lossless or "relaxed" pressure would decrease linearly 
(from p Q ) with x: 

PR = Po - * (36) 


where the overbar indicates a reduced pressure 

P = P/Po u o 2 (37) 

The reduced lossy pressure variation can be gotten by 
integrating equations (33) and (34), with Ap(0) = 0, to get: 


Ap 


T 

T+l 


X + 



T ' 
T+l 


(1+X)-VT 


1 


(38) 
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which yields the interesting relation that volume viscosity will 
cause Ap to be negative, since (l+x) -1 / T - 1 is negative for 
x>0? whereas the pressure lag contributes an additional effect 
which is positive rather than negative. Since Ap = -#e/ R e , when 
T=0, the pressure lag effect changes the lossy pressure behavior 
fundamentally, by eliminating the instantaneous and constant 
Ap contribution, and introducing a contribution which is 
initially zero and becomes positive. 


The significant relaxing degrees of freedom in air are the 
vibration of N 2 , vibration of 0 2 , rotation of N 2 and 0 2 taken 
together, and translation of all constituents. For relatively 
long time scales (>1 ns) only the two vibrational relaxations 
need be considered. Applying equation (38) together with (24a) 
and designating N 2 and 0 2 with subscripts "N" and ”X," we find 
for the reduced pressure: 


P(x) = Pq + X N 


+ X X 


t n _ 

V -f- 


t n 

X T 

t n +i 

R e 

T N +i 

T X - 

. _ v 

'*X 

T X 

X T 

T X +1 

Re 

T X +1 


(1+x) _1 /T n 

(l+x) -1 /Tx 


■]) 

■]} 


(39) 


Values for the relaxation times and relaxation strengths for 
N 2 and 0 2 are taken from reference 6. Then equations (12a, b) and 
(23a, b) are used to find Tp S , r vs , * 7 p, and i/ v . These are listed 
for a temperature of 20°C and humidity values of 0 and 1 mole 
percent in table 1. Because of the small relaxation strengths of 
N 2 and 0 2 , both Tp S and r vs are nearly equal to the 
experimentally determined relaxation time r for each gas 
constituent. 


Equation (39) is plotted in figures 3 and 4 for a temperature 
T = 20°C, length scale L = lm, and absolute humidities h = 0 and 
1 mole percent, respectively. The terms T/(T+1) related to 
pressure relaxation lead to a positive pressure increment, while 
the terms /c/Re related to the traditional volume viscosity lead 
to a negative pressure increment (or decrement) . A net pressure 
decrement is possible only if the volume viscosity terms exceed 
the pressure relaxation terms, as occurs at subsonic Mach 
numbers. In fact, at a speed Mach 0.9 the decrement turns into 
an increment at a reduced distance x * 0.4 - 0.5. At Mach 1 and 
higher the increment is positive over the whole range of reduced 
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distance. Unlike the acoustical case, there do not appear to be 
readily discernible conditions for which the increment reaches a 
maximum . 

Conclusions 

The traditional volume viscosity cannot be given the same 
fundamental status as dynamic viscosity. Not only does the 
traditional volume viscosity fail to evolve naturally from kinetic 
theory , but we have presented explanations showing that the volume 
viscosity coefficient must represent multiple molecular 
relaxation processes, including volumetric dissipation, 
simultaneously when the conventional Newtonian-fluid model is 
used. Furthermore, because of the direct relationship between 
the divergence of the velocity vector and the particle rate of 
change of density, volume viscosity in effect forces normal 
stresses, including pressure, to be modeled in a manner which is 
inconsistent with the more rigorous requirements resulting from 
nonequilibrium acoustic equations of state. We have proposed 
that a pressure relaxation contribution be added directly to the 
constitutive model for normal stresses, in order to bring the 
constitutive model into better agreement with accepted 
nonequilibrium behavior. 

Using accepted acoustical theory, we have estimated the 
values for the pressure relaxation coefficients and volume 
viscosities in pure nitrogen and oxygen at 20°C. Those 
estimates indicate that the nitrogen contributions are nearly an 
order of magnitude larger than the oxygen contributions, but both 
coefficients are influenced strongly by humidity. Our estimates 
show that the volume viscosity is several million times larger 
than the dynamic viscosity for perfectly dry nitrogen, but the 
ratio of volume to dynamic viscosity decreases by a factor of 32 
at moderately high humidity levels (at 20°C) . In addition, if 
pressure fluctuations are on the order of 1 Pa/s, the pressure 
relaxation contribution to normal stresses is nearly an order of 
magnitude larger than the volume viscosity contribution, when the 
modified constitutive model is employed. 

The application to a simple problem in one— dimensional 
convective flow reveals that there are conditions for which the 
volumetric losses in air are negligible (low Mach number, 
moderate humidity), in agreement with the Stokes' hypothesis, but 
other conditions for which they are substantial (high Mach 
number, low humidity). The latter case suggests caution to the 


14 


advocates of the popular notion that computers will replace wind 
tunnels, for a computational algorithm is effective only to the 
extent that it incorporates the significant physical 
interactions. A computation ignoring or misrepresenting 
volumetric losses may fail to predict effects which would be 
observed in wind tunnel tests. 

The constitutive model proposed here should apply equally 
well to liquids, since the acoustical equation of state is the 
same as for gases [7]. The constitutive coefficients can be 
determined from acoustical relaxation data through equations 
(23a-b) . Then the volumetric losses due to thermal relaxation 
[8], for example, can be computed for convective flow in liquids. 
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Table 1. Volume Viscosity Parameters for N 2 and 0 2 at 20°C 


Humidity, mole percent 
Air constituent 
Parameter 

Unit 

0 

n 2 

°2 

1 

n 2 

°2 

Relaxation strength 

€ 


0.00016 

0.0032 

0.00016 

0.0032 

Relaxation times 

T 

/iS 

17700 

6630 

551 

5.37 


T ps 

US 

17700 

6640 

551 

5.38 


T vs 

VS 

17700 

6620 

551 

5.36 

Volume viscosity 

1 p 

/is 

-17700 

-6620 

-551 

-5.36 

coefficients 


Pa. s 

-2512.6 

-943.6 

-78.24 

-0.764 

Mole fraction 

X 


0.79 

0.21 

0.79 

0.21 
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Figure 1. Sound speed ratio squared (c/cq) 2 , showing the Figure 2. Sound absorption per unit wavelength, showing the 

dispersion characteristic of a single relaxation process. relaxation peak characteristic of a single relaxation process. 



Without volume viscosity 
With volume viscosity 


O 
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Figure 3. Reduced pressure versus reduced distance for a linearly accelerating flow in air at various downstream 
Mach numbers. The temperature is 20°C, pressure 1 atm, humidity 0 percent mole ratio, and scale length L= lm. 



Without volume viscosity 
With volume viscosity 
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Figure 4. Same as figure 3 except for a humidity of 1 percent mole ratio (43.4 percent relative humidity). 
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